Time reversal of light by linear dispersive filtering 

^2 ' near atomic resonance 
o 

(N . 

H Linget, T Chaneliere, J-L Le Gouet and A Louchet-Chauvet 

d , Laboratoire Aime Cotton, CNRS UPR3321, Univ. Paris Sud, batiment 505, campus 

univcrsitairc, 91405 Orsay, France 

, E-mail: anne.chauvet@u-psud.fr 

Abstract. Based on the similarity of paraxial diffraction and dispersion mathemat- 
^ . ical descriptions, the temporal imaging of optical pulses combines linear dispersive 



filters and quadratic phase modulations operating as time lenses. We consider pro- 
gramming a dispersive filter near atomic resonance in rare earth ion doped crystals, 
which leads to unprecedented high values of dispersive power. This filter is used in 
an approximate imaging scheme, combining a single time lens and a single dispersive 
section and operating as a time reversing device, with potential applications in radio- 
J>-^ frequency signal processing. This scheme is closely related to three-pulse photon echo 

i-C ' with chirped pulses but the connection with temporal imaging and dispersive filtering 

emphasizes new features. 
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^f) ■ 1. Introduction 

cn 

^ ! The time reversal of a light pulse obeys a very general condition imposed by causality. 

Namely, in order to be ultimately time-reversed, the entire light pulse has to be shelved 
within the time reversing device at some moment. Indeed, no radiative emission should 
^ ! be expected until the pulse tail, to be converted into the time-reversed pulse head, 
has been captured. One cannot help noticing some similarities with other actively 
investigated problems, such as slow light and optical memories. 

In addition to the scientific implications suggested by this general context, the time 
reversal of light also offers stimulating technological prospects. As initially demonstrated 
in the acoustic wave domain, space phase conjugation has to be combined with time 
reversal to concentrate pulsed energy. This has been applied to medical lithotripsy 
techniques pQ. In the framework of electronic warfare one is faced with similar challenges 
in the radio-frequency range. If time reversal is conveniently achieved at low frequencies 
by electronic sampling, recording and playing back [2], the broadband regime requires 
alternative approaches, such as all-optical processing of an upconverted radio-frequency 
waveform. 
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In the present paper we analyze time reversal in the framework of Fourier optics, 
taking advantage of the similarity of linear group delay dispersion and paraxial 
diffraction. Within this picture, we disclose an approximate temporal imaging scheme 
that time reverses the incoming waveform in the same way as a lens generates an inverted 
image. This scheme is considered for time-reversal of optically-carried radio-frequency 
signals. 

As pointed out above, the waveform to be time-reversed has to be entirely shelved 
in the device at some moment. The condition is handled quite easily with pulses shorter 
than a few nanoseconds. The task is more challenging when the waveform duration 
reaches a few microseconds. To manage this issue, we use a linear dispersive filter 
working in the vicinity of atomic resonance. The waveform is saved in long-lived atomic 
superposition states, just enough time until the signal tail is captured and launched 
back as the head of the reversed light pulse. 

In Section [21 we place time reversal within the general framework of temporal 
imaging. After a rapid historical survey we consider an approximate imaging scheme 
that gives rather simple access to time reversal. We point out the connection with the 
ancient pinhole camera. In Section [3] we consider the devising of the linear dispersive 
filters required for temporal imaging. We show that such filtering can be provided 
by appropriately programmed media near atomic resonances. Experimental results are 
presented in Section HJ 

2. Temporal imaging and time reversal 

2.1. Fourier optics and temporal imaging 

With the advent of microwave and optical coherent sources, Fourier optics was actively 
developed in the sixties [3]. Quite soon it was noticed that the paraxial diffraction 
of a spatial wavefront obeys the same mathematical description as the propagation 
of a time domain waveform through a dispersive medium. Both descriptions rely on 
Fresnel transforms that just differ by the transform variable. The transverse space 
coordinate of diffraction is replaced by time in dispersion. This time-space duality was 
first pointed out in the framework of RADAR pulse compression [U [5] . The extension 
to optical domain [6] was related to the rapid progress of laser pulse compression in the 
femtosecond range (TJ, [HJ [9J [10] . 

Typically, pulse compression involves two steps. First, additional spectral 
components are generated by a non-linear process such as self-phase modulation. Ideally 
this is equivalent to "chirping" the pulse frequency, i.e. modulating the incoming field 
with the quadratic phase factor e ir * 2//2 , where r stands for the chirp rate. While 
preserving the temporal profile, the transformation broadens the spectrum, displaying 
the different spectral components in a time-ordered sequence. In a second step, the 
time separated spectral terms are synchronized by linear filtering through a matched 
dispersive element, which results in pulse compression. The dispersive filter multiplies 
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the spectral component at u by the quadratic phase factor expi0(u;) where (j)(cu) = 
HU 2 /2. The dispersion stage is matched to the chirping step if fir = 1. Then the 
emerging waveform reads as £(rt)e irt2 ^ 2 , where £{oj) represents the time-to-frequency 
Fourier transform of the incoming pulse £ (t) . Hence the waveform spectrum is displayed 
in the time domain at the filter output. Compression can occur provided 1/r is smaller 
than the squared duration of £(t). 

The frequency chirping step generates spectral components in the same way as a 
conventional lens adds wavevectors to the spatial wavefront. The time-domain phase 
shift rt 2 /2 is reminiscent to tt(x 2 + y 2 )/Xf, the phase shift that affects a wavefront at 
wavelength A, along the path at transverse position (x, y) through a conventional thin 
lens with focal distance /. Then, the dispersive filtering step with \ir — 1 is equivalent 
to the free space diffraction of a wavefront from the lens to the focal plane. The two-step 
process performs the Fourier transform of the incoming waveform in the same way as 
a lens takes the far field or Fraunhofer diffraction pattern of an incident wavefront into 
the focal plane [TT] . 

Therefore, time space duality in optics is not limited to the similarity of diffraction 
and dispersion. It also extends to the lensing function. Actually the mathematical 
analysis of spatial imaging can be applied to time-domain waveforms, and any coherent 
optics element can be assigned a counterpart in the time domain. From this emerges 
the concept of temporal imaging [T2l [T3] . 

Temporal imaging involves three steps at least. The incoming waveform first 
passes through a dispersive element, with dispersion coefficient fj,i, that operates as 
a region of free-space Fresnel diffraction in the space domain. Then the waveform is 
frequency chirped by factor e* r * 2//2 , which is similar to the action of a thin lens on the 
transverse spatial profile of a beam. Finally, the waveform passes through a second 
dispersive region, with dispersion coefficient /i 2 , equivalent to the lens-to-image free- 
space propagation in space-domain optics. Provided dispersion in the input and output 
sections is correctly matched to the degree of phase modulation by the time lens, a 
time- reversed replica of the input waveform is received at the output. The matching 
condition, 1 / fii+1/ fi2 = r, is similar to the thin lens conjugation formula. The stretching 
factor — /i 2 //ii, depending on the relative strength of the two dispersive sections, replaces 
the space domain magnification factor [T3j. Just as in the space domain, a single lens 
correctly images the object intensity distribution, but alters the phase. A two-lens 
telescope [15] is needed to properly image the field distribution. 

Although most effort has been concentrated for a long time on exactly matched 
imaging or Fourier transform schemes, a new interest has arisen quite recently in 
situations where the Fraunhofer condition is only approximately satisfied. The 
approximate Fourier transform capabilities of a simple dispersion line was first identified 
in the context of real-time spectral analysis [T6l [TT] . It was pointed out that, although 
no lens was used to bring the far field zone to finite distance, the Fraunhofer condition 
could be satisfied provided the squared duration of the incoming pulse is smaller than 
the line dispersion coefficient /x. This procedure, known as wavelength-to-time mapping, 
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was then extended to arbitrary waveform generation [TSJ [B5] and to approximate 
imaging [20|l2T]. 

2.2. Approximate imaging for time reversal 

In many temporal imaging demonstrations, the lens effect is obtained by mixing the 
beam carrying the waveform with a frequency chirped pulse in a two-wave or four-wave 
mixing non-linear material [HI [22} [231 H5]- Without resorting to parametric optical 
processes, one can implement the lens effect in a simpler way by just modulating the 
amplitude of a frequency chirped optical carrier. The price to pay is that the optically 
carried waveform is created by the lens itself. In other words, the optical signal does not 
exist before shaping by the lens. This rules out the input step of propagation through a 
dispersive region, required in the temporal imaging scheme (see section l2TTj) . Modulating 
the amplitude of a frequency chirped carrier is equivalent to having a collimated beam 
pass through an object attached to a lens. Combining this input step with a single 
dispersive section, one can build an approximate time-reversed image of the incoming 
waveform. 




-L L 



Figure 1. Conventional optics image reversal schemes. In upper box (a), the object 
at z = —L is illuminated by a collimated beam. The size of the object smallest detail is 
represented by a. The lens, placed side-by-side with the object, concentrates the energy 
on a AL/er-wide bottleneck in the focal plane (z = 0). The imaging spatial resolution 
on the screen at z = L, of order VXL, is limited by the Fraunhofer condition. The 
lower box (b) represents the pinhole camera. The black dashed lines depict the camera 
impulse response. They represent geometrical rays that are traced from a single point 
on the object to the screen, through the a-wide slit. The shaded zones with different 
shade-levels are intended to show that a small part of the incoming energy is allowed 
by the slit to the screen. 



Before considering operation in the time domain, let us first address the 
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conventional optics counterpart, as sketched in figure QJa). An A-wide object E(x) 
is positioned at z — —L against a converging lens with focal distance L. Let a ^ A 
represent the smallest characteristic distance of variation of the object transmission. 
After crossing the lens and the object, an incoming collimated beam travels to a screen 
located at z = L. A detailed description of the propagation of a spherical wave after 
passing through an object can be found in [23]. We choose a different but complementary 
point of view and analyze the field propagation from z = —L to z = L in two segments. 

In the first segment, from z = —L to z = 0, the wavefront emerging from the lens is 
converted into the far field or Fraunhofer diffraction pattern, brought to finite distance, 
at z — 0, by the lens. The focal spot, of size XL/ a, is nothing but the Fourier transform 
E(u/ (XL)) of the incoming signal E(x), multiplied by the phase factor e "™ 2 /( AL ). 

Traveling beyond the focal plane, the wavefront enters the second segment where 
free space propagation again operates as a Fourier transformer in the Fraunhofer limit. 
From z = to z = L, the Fraunhofer condition is satisfied provided the squared spot 
size (XL /a) 2 at z = is smaller than XL, which leads to XL < a 2 . After undergoing 
two successive Fourier transforms, the intensity transverse distribution at z — L, equals 
\E(—x')\ 2 . According to the Fraunhofer condition and to the a < A assumption, the 
AL/cx-wide focal spot behaves as a narrow bottleneck, much smaller than the object 
size, through which all the incoming energy is channelized. The XL < a 2 condition sets 
the distortion- less imaging capabilities of the device. 

Back to the time domain, let the incoming signal S(t) be injected in the dispersive 
medium after being mixed with the optical carrier by the time lens. When the dispersive 
line coefficient, /i, equals the lens inverse chirp rate 1/r, the signal in the "focal plane" 
of the time lens is nothing but the time-to- frequency Fourier transform S(u) of S(t), in 
accordance with the above conventional optics discussion. Specifically, the " focal plane" 
signal Sfoc(t) reads as S(rt)e trt2 ^ 2 . 

If fi > 1/r, the waveform travels beyond the "focal plane". Just as free space 
propagation in conventional optics, the dispersive line section downstream from the focal 
plane, operates again as a Fourier transformer in the Fraunhofer limit [161 El EE] • When 
fj, = 2/r, the dispersion coefficient of this section reads 1/r, and the Fraunhofer condition 
is satisfied if the squared duration of Sf oc (t) = S(rt)e"' t2 ^ 2 is smaller than 1/r. The 
duration of Sf oc (t) equals (ra)" 1 , where a represents the smallest characteristic variation 
time interval of S{t). Hence the Fraunhofer condition reads as ra 2 > 1 and e irt2 / 2 1 
over the width of S{ oc (t). Finally, after undergoing a pair of Fourier transforms, the 
time-reversed waveform S(—t) is displayed at the output [20] (see figure [2]). 

The image sharpness condition can also be given a group delay picture. Let ujq 
be an arbitrary frequency within the dispersive filter frequency bandwidth. At a given 
frequency u, input and output are separated by the group delay t 9 (uj) = d u (j)(u)). At 
the line output, where /i — 2/r, r g (uj) = r g (uJo) — 2(oj — Uo)/r. Let a Fourier-transform- 
limited, r-long, temporal substructure be injected at frequency u in the dispersive 
medium. This pulse spreads over a spectral interval of order 1/r. Due to group delay 
dispersion d u T g (oj) = —2/r, the pulse undergoes stretching of order \d w T g (txi)\/T = 2/(rr) 
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Figure 2. Approximate imaging for time reversal. Passing through a time lens, the 
incoming S(t) waveform is changed into 5 M=0 (i) = S{t)e irt2 / 2 . Then S^t) propagates 
through a dispersive filter, from [i = at the input to /x = 2/r at the output. The time 
variation of S^(t) is displayed in 3D representation at different //-positions, illustrating 
the waveform compression into the focal plane, at fj, = 1/r, the subsequent stretching 
and the time-reversal at fi — 2/r. The phase evolution in the focal plane results 
from the time separation of the two incoming peaks. The output profile represents 
|<S/i=2/r(0|- I n this simulation, r has been set equal to 20/r 2 , where r represents the 
full width at half maximum of the incoming peaks. 



while travelling through the dispersive filter, and preserves its initial duration provided 
t > 2/(rr). Hence, J^/r represents the duration of the shortest substructure that 
can propagate throughout the medium without distortion. The proper time-reversal 
bandwidth is thus limited to Jr/2. 




Although attractive in the time domain, this scheme is never considered in 
conventional optics. Indeed, as soon as a lens is available, one takes care of conjugating 
the object and the image through the lens, which relaxes the distortion-less bandwidth 
condition. However, the depicted approximate imaging design is rather reminiscent to 
the pinhole camera, an ancient lens-less imaging device. 

2.3. Connection with the pinhole camera 

To draw the parallel with the pinhole camera, we return to conventional, space domain, 
optics. In the above description we pointed out that the focal spot can be regarded as 
a narrow bottleneck through which all the wavefront information is conveyed. 

Channelizing through a narrow opening obviously evokes the pinhole camera - a 
light-proof box with a small hole on one side. Light from a scene enters through the tight 
opening and projects an inverted image on the opposite side of the box [see figure [11(b)]. 
The pinhole camera operates in Fraunhofer conditions, just as the device shown in 
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figure QJa), and offers the same distortion-less resolution. To satisfy the Fraunhofer 
condition, the opening should be smaller than « V XL. To minimize the low-pass 
filtering of angular frequencies, one should not reduce the opening size beneath this 
value that appears to be the best trade-off between the Fraunhofer and the diffraction 
limits. 

Therefore the two devices offer the same distortion-less resolution, as imposed by 
the Fraunhofer condition. However, despite their similar geometry, they do not share 
any additional common property. 

First, the Fraunhofer condition is always satisfied inside the pinhole box, as soon 
as the opening is properly sized. When object details are smaller than vAL, they 
give rise to high angular frequency components that are just cut-off by the hole. In 
the image on the back side of the box, the details are simply blurred by this low-pass 
filtering effect. In the figure HJa) device, the high angular frequencies are not cut but 
they violate the Fraunhofer condition. We shall examine this point in depth further on. 

Second, the pinhole camera poorly operates in coherent illumination conditions. 
Indeed, if a a-wide detail on the object, located at distance x from the hole axis, is 
illuminated by a collimated coherent beam, the transmitted light travels in forward 
direction within a diffraction limited aperture. Hence it reaches the hole only if the 
~ AL/er-wide projection is smaller than x. Otherwise, the rays are just stopped by 
the light-proof wall of the box. Hence, only a small fraction of the object, very close 
to the hole-axis, can contribute to the image on the back side of the box. This issue 
can be solved by spatially incoherent illumination, with a correlation distance along Ox 
much smaller than a. This sub-modulation gives rise to the angular frequencies that 
are needed to have some light reach the hole from any point of the scene. However, the 
incoherent illumination of the scene generates angular frequencies in a random way, so 
that most of the incoming light is lost, as illustrated in figure QJb). On the contrary, 
the figure []Ja)-device correctly works with coherent light since the lens generates the 
proper angular frequencies to direct all the light through the bottleneck. 

In a time domain pinhole camera, the waveform to be time-reversed would modulate 
a broadband carrier, with an autocorrelation time much smaller than the waveform 
shortest temporal detail. A short temporal gate would play the part of the pinhole, 
through which dispersed components of the initial waveform would be transmitted [25] . 

3. Time reversal through a shaped-absorption medium 

3.1. Linear dispersive filtering near atomic resonance 

Any transparent medium can be used as a dispersive filter. Since the dispersive effect 
grows with the material thickness, hundred-of-meter-long optical fibers can be regarded 
as emblematic media for such applications. Specific dispersive properties are usually 
intrinsic to a given material. Here we consider a different class of dispersive filters that 
are programmed in an absorbing material. 
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Figure 3. Spectral profile of the absorption coefficient as defined by equation( [T]) 
(upper box) and of the corresponding index of refraction variation 5n(uS) = ra(aj)— n(uj$) 
(lower box). When y/JI(io — loo) is large enough, a(u) and 5n(uj) exhibit the same quasi- 
sinusoidal shape with the same period and oscillate in quadrature. 



Let the absorption coefficient of a linear filter be given the following spectral shape: 

a(u) = a (1 + (3 cos [<j>(u - w )]) /(I + P) (1) 

where < /? < 1, Uq is an arbitrary frequency origin, and: 

= fiu 2 /2 (2) 

with yU > 0. According to equation (JTJ, a(u) oscillates non-periodically between a an d 
«o(l — + P), as illustrated in figure El The spectral shape is uniformly engraved 

throughout the depth of the filtering medium. Hence a(u) does not depend on the 
depth coordinate z. 

Let the field £{z = 0,u) be incoming into the filtering material at position z = 0. 
At arbitrary depth z, the transmission factor reads as: 

J~^ z w ) _ e -ikn(uj)z-±a(uj)z _ Q -ik [n(w )+3X(w)] * ^ 

where n(u) and xi^) respectively represent the index of refraction and the electric 
susceptibility Those quantities are linked together by the Kramers- Kronig relations 
that express causality. Term-by-term identification of the exponents in equation ([3]) 
leads to: 

X(u) = 2 [n(w) - n(w )] - ~a(u)) (4) 
According to the Kramers-Kronig relations: 

n(u) - n(u ) = -^-H{a}(u) (5) 
where the Hilbert transform H{f}(x) is defined as: 

H{f}(x) = -P r ^dx' (6) 

71 J — co X — X 
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One easily shows that, provided 4>(u — u ) ^> 1: 

H { e ±^(<"-<"o)} (w) = t* e^^-^, (7) 



which leads to: 

i «o 



1 + (3 e -^~"o) ( 8 ) 
k 1 + L J v 7 

Substituting equation (El) into equation (j3J) and expanding the exponential in series one 
obtains: 



where 



T{z,lo) = %n(z,u) (9) 

m=0 

TmM = i M^y e -f^-^--)-^ (10) 

m! \ 2 1 + p y 

The m = component of the series accounts for direct transmission. The next 
terms behave as the transmission factors of dispersive lines. Each dispersive line is 
characterized by the group delay: 

T^ m \uj) = md(j)(oj — uo)/du = mfi(uj — ljq) (11) 

that grows linearly with u. In the following, for the sake of simplicity, t^(lo) = T g (co). 

Let us concentrate on the first dispersive term, with m = 1. The inverse group delay 
is nothing but the period of the absorption coefficient spectral modulation. Hence, the 
maximum size of the group delay is determined by the narrowest spectral feature one 
can engrave in the filtering medium. 

Affected by different group delays, the successive dispersive terms can give rise to 
temporally separated contributions. This will be made clear in the following. If the 
successive contributions do not overlap in the time domain, the transmission efficiency 
of the mth term can be defined as: 

r] m (a z,(3) = \T m (z,uj)\ 2 . (12) 

One easily verifies that: 



ElU<**fl =•-»!. (£*f) (13) 

where Io(x) stands for the 0-order modified Bessel function. The m = 1 response 
efficiency reads as: 

"(sr?!)' 6 "* (14) 

The maximum efficiency for the first-order response, 13.5%, is reached at a^z = 4, f3 = 1. 
The same fraction of the incoming energy is carried away by the direct transmission 
component: r)o(aoz,(3) = e _T +?. According to equation ffl3|) . all the higher order terms 
collect less than 4% of the incoming energy. Delayed transmission through a spectrally 
periodic sinusoidal filter is commonly known to be limited to 13.5% efficiency. Not 
surprisingly, the locally periodic approximation leads here to the same figure. 
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Dispersive filtering near atomic resonance has already been used in combination 
with a time lens for real-time spectral analysis [26 1 |2T | 128] and analog arbitrary waveform 
generation [291 EQ] in the prospect of RADAR applications. 

When addressed in the context of atomic resonance, both the absorbing material 
programming and the signal filtering are usually identified as parts of a stimulated 
photon echo process (31]. Photon echo is known as a four wave mixing process, and 
it may be misleading to refer to non-linear optics when one really deals with linear 
filtering. The first two pulses of the photon echo sequence are just used to program 
the filter. Reference to photon echo might also restrict the filter programming to the 
canonical two pulse sequence. As a matter of fact, a vast variety of techniques has 
been used to program the material, as illustrated by various demonstrations ranging 
from the accumulated photon echo [32J, [33] to the rainbow analyzer [34] and the atomic 
frequency comb (AFC) [351 EEl E7]. Because of the flexibility of the programming step, 
we purposely let aside the specific preparation technique, concentrating on the most 
general features that can give us some insight on ultimate properties and limitations. 



3.2. Application to time reversal 

In this section we derive the image of an incoming waveform through our approximate 
imaging scheme in a single step, following the calculation achieved by J. Azana et al. in 
the spatial domain [23] and in the time domain [2"T] . 

Let the waveform S(t), centered at t — 0, be upconverted on the frequency chirped 
optical carrier C(t), resulting in an input field £{z = 0,t) = S(t)C(t). The carrier field 
reads as: 

C(t) = C e iuJlt ^ rt2 (15) 

The instantaneous carrier frequency uj\ — rt equals uji at t = and decreases with time. 
Within the waveform duration T, the carrier frequency satisfies the dispersive filter 
operation condition <f>(u — oj ) ^> 1, which leads to: 

rT/2 < u)i - u (16) 

Propagating through the m — 1, L-thick dispersive filter, the input field £(z = 0,t) 
undergoes the transformation: 

£(L, u) = Ke-^ iuJ - uo)2 - lkL £{0, u) (17) 



where k = yr)i(otoL, f3), and where the time-to-frequency Fourier conjugation is defined 
by: 

/» = J dtf(t)e~^\ fit) = i- J dwf(u)&*. (18) 

Let r/i = 2, which is the matching condition for time reversal with unit 
magnification. A tedious but straightforward calculation leads to: 

£(L,t) = KC e im ^- [ dooSi^e-^'^ 1 ^^ 2 (19) 
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where: 



tf>{t) = u ± t + -r [t - Tg^)} 2 - -rr 9 (wi) 2 



(20) 



If, in the integrand of equation f|T9|) . fiu 2 <C 2tt over the waveform bandwidth, the sum 
over u reduces to the inverse Fourier transform of S(u>) at time —t + T g (ui) and the 
transmitted field can be expressed as: 

£{L, t) = kCq^S (-t + r fl (wi)) (21) 

As pointed out above, the ignored phase contribution, fiu 2 , reflects the imaging 
imperfection of the approximate imaging scheme. 
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Figure 4. Carrier and signal time diagrams. The incoming carrier instantaneous 
frequency decreases with time from cj max to w m j n . On the contrary, the first-order 
response carrier frequency increases with time from w m in to u; max . The (w m in, Wmax) 
interval is centered at u)\. The outgoing response, centered at t = T g (oj\), is time- 
reversed with respect to the incoming signal, centered at t = 0. The carrier frequency 
range [w m ; n ,cj max ] and the signal duration T are connected by 2T = T g (oj max ) — 
T g (w m in). In the figure, the first (m = 1) and second (to = 2) order responses partly 
overlap. To avoid this, the signal duration must satisfy the condition T < r g (w{)/2. 



The time-reversed waveform S (— t + r g (oJx)) is carried by the field Coe 1 ^^ . 
According to Eqs. ( !T5|) and ( )20|) . the sign of the carrier chirp rate at the output is 
reversed with respect to the input (see figure BJ. 

Finally the restored waveform should not be spoilt by higher order responses, as 
noticed in section [37T1 As illustrated in figure (jlj), the first order response ends at time 
tmox = r g( w i) + T/2. The second order response starts at t^ in = t^\uji) — 3T/2 = 
2r g (ui) — 3T/2. Those responses do not overlap temporally, provided t$ ax < t^ in , 
which reduces to: T < T g (ui)/2. The carrier frequency range [w m i n , o; max ] and the signal 
duration are connected by 2T = r ff (u; max ) — r g (u; m i n ). 
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4. Experimental 

In this section we present the experimental implementation of time-reversal with the 
approximate scheme described in Sec I2.2[ where the dispersive filter is created in an 
absorbing medium. 



4-1. Preparation of a dispersive filter 

Famous for its absorption line in the telecom wavelength range, the erbium Er 3+ 
ion inserted in yttrium orthosilicate (Y 2 Si05) offers spectroscopic properties allowing 
preparation of a dispersive filter with a large dispersive power. On the one hand, each 
individual Er ion has a homogeneous linewidth of a few kHz when cooled down at 2 K 
in a 2-tesla magnetic field. On the other hand, environment inhomogeneity at the ion 
sites induces an inhomogeneous distribution of the individual lines, resulting in a GHz- 
wide absorption profile for the ensemble. Shining a spectrally narrow laser on such a 
profile transfers resonant ions to their excited state and creates a long-lived dip in the 
absorption profile. This spectral hole-burning process is at the heart of the preparation 
of our dispersive filter. 

The preparation of the filter is inspired from the usual stimulated photon echo 
scheme. The first two pulses £i(t) and E^it) of a stimulated photon echo sequence burn 
a sinusoidal spectral grating in the absorption profile, with a spacing equal to the inverse 
time separation of the two pulses (figured]). One can extend this observation to a pair 
of chirped pulses with opposite chirp rates — rp and rp: each homogeneous frequency 
class of atoms senses a pair of pulses with a frequency-specific time separation [271 12H] • 
This way, the inhomogeneously broadened profile is shaped into a quasi-sinusoid with 
linearly varying period, hence the dispersive filter with \i = 2/rp (figure [6]). 




Figure 5. Two monochromatic pulses separated by t\i induce a sinusoidal modulation 
of the absorption profile with spacing 2ir/ti2- 



The crystal used in our experiments is a 10 mm-long 0.005% Er:YSO crystal, cooled 
down to 2 K in a helium bath cryostat. A 2-tesla magnetic field is applied by means of 
superconducting split coils in the sample space. 

We prepare the dispersive filter with two oppositely chirped pulses sent to the 
Er:YSO crystal, sweeping their frequency over A/2tt = 1.09 GHz in 6 /xs, resulting in a 
chirp rate rp = 1.14 x 10 15 rad.s -2 . We probe the modified absorption profile by shining 
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Figure 6. Two chirped pulses with opposite chirp rates — rp and rp induce a quasi- 
sinusoidal absorption grating with a variable spacing. 



an attenuated, slowly chirped laser pulse (2 MHz in 40 /is) around three frequencies (see 
figure [TJfor details). We use a single laser source for the preparation pulses and the probe 
pulse, namely an agile extended cavity diode laser that has been designed to perform 
reproducible GHz- wide frequency sweeps within a few /xs without mode hops [SB]- The 
laser frequency sweeps are obtained by feeding a linear voltage ramp to an intra-cavity 
electro-optic crystal, with smooth slope changes to avoid excitation of piezo-electric 
resonances in the crystal. The laser is fed into an erbium-doped fiber amplifier and 
shaped in amplitude with a free-space acousto-optical modulator. 

The transmitted probe intensity is shown in figure [71 As expected, we observe 
a locally sinusoidal modulation with frequency- dependent spacing, corresponding to 
the pulse separation sensed by the atoms absorbing at that frequency. The resulting 
dispersive power of such a filter is/i = 2/rp = 1.5 ms/nm, that is 8 orders of magnitude 
larger than the dispersive power of a standard km-long fiber at 1.55 /xm (17 ps/nm/km). 
The measured contrasts /3 range from 0.08 to 0.15. 



(9 = 0.15 (3 = 0.10 /3 = 0.08 




Probe frequency detuning 5 (MHz) 



Figure 7. Normalized transmission measurement around different frequencies within 
the 1.09 GHz engraved frequency range (blue line). After the two programming pulses, 
a weak slowly chirped pulse is sent to the crystal and its transmitted intensity is 
monitored. Depending on the spectral range probed, one observes a locally sinusoidal 
modulation of the transmission, with frequency-dependent spectral period l/t\2 and 
contrast /3. Thin black line: [1 + /3 sm(2-K5t r2 + 6)]/(l + (5). 



We insist on the fact that the 13.5 % maximum efficiency announced in section [37X1 
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corresponds to a quasi-sinusoidal grating with unity contrast and a maximum optical 
depth OqL = 4. According to equation [HI the efficiency in the first order response varies 
like (3 2 . With an opacity a^L = 2 and based on the contrast measurements shown in 
figured we expect a retrieval efficiency 771 ~ 1.3%o. 

The low contrast achieved in our experiments originates from the perturbative 
character of the preparation step. However, thanks to its short duration (~ 15 /xs), 
the preparation step can be repeated many times within the grating lifetime, similarly 
to accumulated photon echoes [32], [33] . This should increase the grating contrast and 
therefore enhance the efficiency, up to 13.5%. Besides, optimization of the grating 
shape is known to allow even higher efficiencies beyond the 13.5% limit. For example, 
an efficiency close to 20% has been demonstrated with a rectangular grating in Tm:YAG, 
in the context of quantum memories [39|. 

4-2. Time-reversal of a 5 fis-long signal 

The complete pulse sequence is shown in figure [HJ To ensure time-reversal with unit 
magnification, the signal S(t) is upconverted on a chirped optical carrier with chirp 
rate r$ = —rp = — r with an acousto-optic modulator. The signal duration (5 fis) 
is compatible with typical radar signal durations. The signal propagates through the 
dispersive filter programmed in the Er:YSO crystal by the preparation pulses. The first 
order response of the filter is emitted as a chirped sequence of pulses with rate r. 




Time (|as) 

Figure 8. Experimental pulse sequence. Thick lines: incoming preparation and signal 
pulses. Thin line: dispersive filter first order response, emitted as a chirp with rate r, 
carrying the time-reversed signal 

We use the experimental setup described in section 14.11 to generate the whole 
sequence. An additional acousto-optical modulator attenuates the strong pulses 
transmitted through the Er:YSO crystal, to protect the avalanche photodiode Thorlabs 
APD110C used for detection. 
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We emphasize that the absorption profile of our Er:YSO crystal has a full width at 
half maximum of 2.7 GHz, and its optical depth a^L is not constant across the 1.09 GHz 
dispersive filter width. The filter contrast f3 is not constant either across the filter width, 
as shown in figure [71 A variation of the efficiency in the m = and m = 1 responses 
should therefore be expected [see equations (flO]) and (TT21)] . 

Given the measured grating contrast (f3 ~ 10 %) and typical optical depth in our 
sample a L = 2, we expect a first-order response efficiency r]i ~ 1.3%o, and a second 
order response much less efficient than the first order one: 772 = 2 x 10~ 3 ?7i. This is why 
we allow a time gap between the signal and the echo that makes the second order partly 
overlap with the first order response (figure Hj). 

Due to causality, the whole input signal must be confined before the time-reversed 
response is emitted. During this time, the information is saved in the crystal as 
a superposition state that is sensitive to decoherence. The various mechanisms at 



play are described in Appendix A With our experimental parameters, the decoherence 
mechanisms remain small enough and the time-reversal efficiency is not affected. 

The time-bandwidth product P of our time-reversal protocol is limited by the 
intrinsic properties of the atomic ensemble. Indeed, P reads as: T x 5 where T is the 
input signal duration and 5 is its spectral width. The Fraunhofer condition imposes that 



P <C yAT/2. A cannot exceed the inhomogeneous linewidth of the atomic ensemble. 
With a signal duration of 5 yus, this leads to ? < 80. Using a different crystalline 
host such as LiNbOa, where the Er inhomogeneous linewidth reaches 180 GHz and the 
decoherence effects are similar to that in YSO [30], would allow a larger bandwidth and 
a larger time-bandwidth product for our time-reversal protocol. This would require a 
laser able to perform fast reproducible sweeps over such a frequency range. However, as 
long as approximate imaging is used, the Fraunhofer condition remains. 

4-3. Exploring the Fraunhofer condition 

In this section we investigate the signal distortion due to the violation of the Fraunhofer 
condition. We explore the effect of the Fraunhofer condition by sending signals with 
different bandwidths on the dispersive filter described above. All the incoming signals 
consist in a train of Gaussian pulses J2i ex P — 2^ with r ranging from 18 to 67 ns. 
Successive peaks are separated by multiples of 207 ns. The Fraunhofer condition is 
satisfied when r 3> y/2/ (2iiy/r) = 17 ns. 

Figure [9] shows the transmitted signals and their first-order response detected on 
the same avalanche photodiode. The transmitted signal is the n = response of the 
dispersive filter. Its intensity varies because of the crystal optical depth a^L and filter 
contrast /3 variation over the filter frequency range, as expected from our analysis of 
figure [71 

The m = 1 response contains the time-reversed counterpart of the input signal. In 
these experiments, the maximum efficiency of the time-reversal process reaches 1.6%o, in 
agreement with the measured population grating contrast. Like the m = response, the 
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Figure 9. Around the Fraunhofer limit. Thick lines: experimental transmitted signal 
and first-order response measured on the same detector after propagation through the 
dispersive filter described in figure [7] The input signal consists in a train of identical 
Gaussian pulses with typical widths r ranging from 18 to 67 ns. The variable intensity 
of the experimental m = and m = 1 responses is due to the variation of optical 
depth and filter contrast across the filter frequency range. Thin lines: input signal 
and calculated response intensity of an ideal dispersive filter with dispersive power 
jjL = i.5 ms/nm, only adjusted in height to be compared easily with the experimental 
signal. The curves are vertically offset for clarity. 



first-order response intensity varies over the signal duration because the optical depth 
OqL and the filter contrast (3 are not constant over the filter frequency range. 

The signal is faithfully reversed as long as r > 40 ns, in agreement with the 
Fraunhofer condition r ^> 17 ns derived above. When the signal spectral width exceeds 
the filter bandwidth, each peak of the time-reversed signal is stretched and overlaps 
with its nearest neighbours that are emitted at different frequencies, hence the beatnote 
between close pulses. As expected, high frequency components in the input signal induce 
strong distortion in the filter response. We compare the experimental measurements 
with the calculated signals as given by equation ( 119]) and see almost perfect matching. 
We visualize the fidelity of the time-reversal by plotting the correlation coefficient 
between the input signal and its experimental and calculated time-reversed replicas 
in figure [101 



4-4- Simulation of a pinhole time camera 

We calculate the image obtained with a time pinhole camera setup [25], with two 
successive dispersive lines with \i = 1/r = 0.75 ms/nm on either side of a rectangular 
time window. We set the pinhole time window to £ = 40 ns in order to make the 
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Figure 10. Fidelity of the time-reversal with respect to the signal spectral width. 
Circles (resp. solid line): correlation coefficient between the input signal and its 
experimental (resp. calculated) time-reversed replica shown in figure [5] Dashed 
line: correlation coefficient between the input signal and the calculated equivalent 
pinhole camera image. The two setups offer the same bandwidth, with a cutoff around 
(271-t)- 1 = 4 MHz, ie t = 40 ns. 



pinhole camera and the time-reversal setups comparable. Note that the Fraunhofer 
condition is satisfied inside the dark chamber. In the conventional optics domain, 
incoherent light is needed at the input of the pinhole camera to allow some light to 
reach the pinhole even when coming from the outer part of the object. In our time 
pinhole camera simulation, we introduce incoherence in the input signal by multiplying 
it with a random phase distribution factor. The calculated outputs are averaged over 
500 random phase distributions, normalized and shown in figure [TTJ One recognizes the 
pinhole outputs as the convolution of the time-reversed input signals with the impulse 
response of the pinhole setup, ie the Fraunhofer diffracted image of the pinhole aperture: 
sine 2 (r£t/2). This impulse response has a full width at half maximum of 125 ns, larger 
than or comparable to the input signal details, which results in a broadening of the 
signal details and a contrast reduction in all the time-reversed images shown in the 
figure. 

In figure [10] we display the fidelity of time-reversal with a pinhole time camera. 
Perfect time-reversal with a pinhole occurs only when the input signal details are larger 
than the impulse response typical temporal width, ie when the signal spectral width 
(2ttt)^ 1 <C 1.3 MHz. Beyond this value, the pinhole camera fidelity decreases at a slower 
pace than our approximate time-reversal scheme. Indeed, beyond the Fraunhofer limit, 
high-frequency components are simply cut-off by the pinhole, unlike in the time-reversal 
experiment where they strongly distort the reversed signal. The two architectures 
nevertheless share a similar bandwidth, close to j2/r. 




5. Conclusion 

To sum up, we have demonstrated the time-reversal of optically carried signals with an 
approximate temporal imaging scheme, involving a time lens and a dispersive filter. The 
dispersive filter is obtained with a rare-earth-ion-doped crystal whose absorption profile 
has been modified with spectral hole burning. The dispersive power of our filter exceeds 
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Figure 11. Simulation of a pinhole time camera. Left curves: input signals (identical 
to those used in our time- reversal experiments). Right curves: calculated pinhole 
response. The curves are vertically offset for clarity. 

that of conventional fibers by many orders of magnitude. Because of its approximate 
character, our time-reversal scheme suffers from a bandwidth limitation, analog to 
the Fraunhofer condition in optics. We investigate this condition experimentally, and 
observe strong signal distortion when the condition is not satisfied. We connect our 
time-reversal scheme with an equivalent time domain pinhole camera. We observe that 
although they share the same bandwidth, they distort the input signals differently. 

For better time reversal bandwidth, taking full advantage of the filter bandwidth, 
one cannot dispense with true temporal imaging, which involves a time lens and two 
dispersive lines, respectively located upstream and downstream from the lens. Those 
elements are needed to conjugate the temporal object and its time-reversed image. To 
operate with two dispersive elements, the time lens must deal with an optical signal 
at the input and the output. A wave- mixing-based non-linear optical device commonly 
satisfies such conditions [HJ |22} |23j [15]. When the dispersive filters are programmed in 
an absorbing medium, one is faced with the additional requirement of having both filters 
work at the same wavelength. Then the time lens shall be based on difference frequency 
generation (DFG) rather than on more usual sum frequency generation (SFG). For 
instance, entering the lens at 1.5 /im, the signal shall emerge at the same wavelength 
if the lens works as a DFG-mixer with a chirped field at 750 nm. Those non-linear 
techniques are well documented |JT]. 

More problematic is passing the signal twice through the shaped-absorption-based 
filter. According to section l3~TTj single pass optimized efficiency does not exceed 13.5%, 
which leads to less than 2% efficiency through two successive dispersive filters. However, 
we know that the 13.5% barrier can be broken. For instance, optimizing a periodic filter 
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in Tm 3+ :YAG, we have been able to demonstrate about 20% efficiency with rectangular- 
shaped spectral grooves in the context of quantum memory for light |39j . 
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Appendix A. Sources of decoherence 

The efficiency of our time-reversal filter is affected by decoherence mechanisms, 
originating from the finite homogeneous linewidth, but also from spectral diffusion and 
instantaneous spectral diffusion |42j . 

Spectral diffusion refers to the fluctuations of each ion transition frequency as a 
result of interaction with the environment. In crystals doped with Kramers ions, such as 
Er 3+ , spectral diffusion is mainly caused by magnetic dipole-dipole interaction between 
the active ions themselves. This contribution can be strongly reduced by an external 
intense magnetic field that locks the spins in a fixed orientation. In Er:YSO we apply 
a 2-tesla field in the (D\, Z^-plane, 135° from the Di-axis |33], while maintaining the 
sample at 2K. The subsisting spectral diffusion affects three pulse photon echoes (3PE) 
in two different ways. 

On the one hand, spectral diffusion destroys the atomic coherence that carries 
optical excitation during the t±2 interval between the first two pulses. Indeed the 
frequency fluctuations alter the relative phase of the atomic states in a cumulative way. 
The same coherence deterioration takes place between the third pulse and the echo. 

On the other hand, spectral diffusion tends to erase the periodic spectral structure 
that is engraved in the absorption profile by the first two pulses. The signal is the more 
affected as the structure period, given by l/ti2, is narrower, and as the waiting time TV 
from the second to the third pulse is longer. However, the decay saturates as a function 
of T\y since the frequency excursion range is limited. 

In Er:YSO, the observed 3PE intensity ^-dependent decay is consistent with 
exp[— 47rti 2 r e //(ti2, T w )] where the effective decay rate T e ff{ti 2 ,T w ) reads as [4~2] : 



In this expression r , T$d and R respectively represent the single ion homogeneous 
width, the frequency fluctuation range and the inverse correlation time of the 
fluctuations. The T^-dependence saturates when RTw > I. In our crystal we have 
measured r = 800 Hz, T SD = 3 kHz, and R = 10 ± 1 x 10 3 s _1 . 




(A.l) 



Time reversal of light by linear dispersive filtering near atomic resonance 



20 



Our time-reversed signal can be analyzed as a combination of 3PE sequences 
spreading over a range of different i 12 and TV values. With a minimum Tjy value 
of ~ 100 /is we operate in the saturation regime, and with tu < 6 /is we get 
exp[-47rt 12 r e// (t 12 ,7V)] > 0.93. 

Instantaneous spectral diffusion refers to random frequency shifts of each ion caused 
by excitation- induced changes in the ions environment. As ions in the lattice are 
excited by the incident light, changes in the excited ions properties alter the resonance 
frequencies of nearby ions. Based on |Hj we estimate that the preparation pulses induce 
a broadening of the single ion homogeneous linewidth T up to 80 kHz, when all the 
erbium ions of the filter frequency range are transferred to their excited state. In our 
experiments, we minimize the contribution of the instantaneous spectral diffusion by 
choosing a large beam waist diameter (130 /im) in the crystal, a moderate beam power 
(5 mW), and a rather short preparation duration (two 6 /xs-long pulses). 
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